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ABSTRACT 


The seasonal sea circulation and thermohaline structure in the Japan/East Sea 
(JES) were studied numerically using the Princeton Ocean Model (POM) with horizontal 
resolution varying from 11.54 to 18.53 km and 15 sigma levels conforming to a relatively 
realistic bottom topography. A twenty four month control run was performed using 
climatological monthly mean wind stresses and heat and salt fluxes as surface forcing and 
observational oceanic inflow/outflow at open boundaries. The seasonally averaged 
effects of isolated forcing terms are presented and analyzed from the following 
experiments; 1) non-linear effects removed, 2) no lateral transport at open boundary, and 
3) wind effects removed. This procedure allowed analysis of spatial and temporal 
contributions of the isolated parameter to the general hydrology of the JES and some of 
its specific features. Major currents are simulated reasonably well compared to 
observations. The nonlinear advection does not affect the general circulation pattern 
evidently, but does affect the formation of the mesoscale eddies, especially the Ulleung 
Basin (UB) eddy (all seasons) and the Japan Basin (JB) cyclonic gyre (spring). The 
lateral boundary forcing enhances (weakens) the JES volume transport in the summer 
(winter). The wind forcing is the most important factor (80 %) for the generation of tiie 
JB cyclonic gyre. It drives the Liman Current and damps the East Korean Warm Current 
in the winter, and generates the UB eddy, and eddies along the Japan Coast Current 
(JCC) in all seasons. However, it has almost no effect on the JCC for all seasons. 
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I. INTRODUCTION 


The Japan Sea, known as the East Sea in Korea, has a steep bottom topography 
(Fig. 1) that makes it a unique semi-enclosed ocean basin overlaid by a pronounced 
monsoon surface wind. The Japan/East Sea, hereafter referred to as JES, covers an area 
of 10^ km^. It has a maximum depth in excess of 3,700 m, and is isolated from open 
oceans except for small (narrow and shallow) straits, which connects the JES to the 
Pacific Ocean. The JES has great scientific interest as a miniature prototype ocean. Its 
basin-wide circulation pattern, boundary currents, polar front, mesoscale eddy activities 
and deep water formation are similar to those in a large ocean. 

The JES physical oceanography has been investigated for several decades. The 
warm Tsushima Current, dominating the surface layer, flows in from the East China Sea 
through the Tsushima/Korean Strait and carries warm water from the south up to 40°N 
where a polar front forms (Seung and Yoon, 1995). The cold Liman Current flows in 
from the Okhotsk Sea through the Tatar Strait and carries cold water from the north. 
Most of the nearly homogeneous water in the deep part of the basin is called the Japan 
Sea Proper Water (Moriyasu, 1972) and is of low temperature and low salinity. Above 
the Proper Water, warm and saline water, which enters through the Tsushima Strait, 
flows northeastward and flows out through the Tsugaru and Soya Strait. 

The seasonal variability of the JES sea surface temperature (SST) has been 

studied by many investigators (Isoda and Saitoh, 1993; Isoda et al., 1991; Isoda, 1994; 

Kano, 1980; Maizuru Mar. Observ., 1997) using limited data sets. For example, Isoda 

and Saitoh (1993) analyzed the satellite infrared (IR) images in the western part of the 
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JES and the routine hydrographic survey completed by the Korea Fisheries Research and 
Development Agency in 1987 and found the sea surface temperature (SST) patterns in 
winter and spring. They found that a small meander of a thermal front originates from 
the Tsushima Strait near the Korean coast and gradually grows into an isolated warm 
eddy with a horizontal scale of 100 km. The warm eddy intrudes slowly northward from 
spring to summer. Recently, Chu et al. (1998a) reported the seasonal occurrence of JES 
eddies from the composite analysis on the U.S. National Centers for Environmental 
Prediction (NCEP) monthly SST fields (1981-1994). For example, they identified a 
warm center appearing in later spring in the East Korean Bay. Chu et aL (1999) further 
reported the seasonal variation of the thermohaline structure and inverted circulation 
from the Navy's unclassified Generalized Digital Environmental Model (GDEM) 
temperature and salinity data on a 0.5° x 0.5° grid. 

The JES circulation is extremely complicated owing to its variable and diverse 
processes such as annual and interannual variations in lateral transport at straits, wind 
stress, branching of the Tsushima Warm Current (TC), western intensification of the East 
Korea Warm Current (EKWC), flow separation and eddy shedding and sub-polar front 
meandering. Present understanding of the circulation in this region is still in the early 
stage. The recent field observation program on Circulation Research of the East Asian 
Marginal Seas (CREAMS) (Takematsu et al., 1994; Takematsu et aL, 1996) provides 
information on current variability and water mass distribution for the bottom and 
intermediate waters. 

Recently, Chu et al. (1999) used the Navy’s Generalized Digital Environmental 

Model (GDEM) climatological temperature and salinity data on a 0.5° x 0.5° grid to 
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investigate the seasonal variabilities of the Japan/East Sea (JES) thermohaline structure 
and circulations. The GDEM for the JES was built on historical (1930-1997) 136,509 
temperature and 52,572 salinity profiles. A three-dimensional estimate of the absolute 
geostrophic velocity field was obtained from the GDEM temperature and salinity fields 
using the P-vector method. The climatological mean and seasonal variabilities of the 
thermohaline structure and the inverted currents such as the Polar Front, the mid-level 
(50-200 m) salty tongue, the Tsushima Current and its bifurcation, are identified. 

Numerical studies on the JES circulations started in early 1980. Various types of 
models were used such as the multi-layer model (Sekine, 1986, 1991; Kawabe, 1982; 
Yoon, 1982a,b; Seung and Nam, 1992; Seung and Kim, 1995), the Modular Ocean Model 
(MOM) (Kim and Yoon, 1994; Holloway et al., 1995; Kim, 1996; and Yoon, 1996; 
Yoshikawa et al., 1999), the Miami Isopycnal Coordinate Model (MICOM) (Seung and 
Kim, 1993) and the Princeton Ocean Model (POM) (Chu et al., 1998b, 1999a-g; Guo, 
1999; Ro, 1999). Most of the numerical efforts are concentrated on simulating basin¬ 
wide circulation, the TC bifurcation, and formation of the intermediate waters. However, 
there is no modeling study on the formation of the seasonal variabilities of the JES 
circulation and thermohaline structure. 

In this study, we use the Princeton Ocean Model (POM) to investigate 

mechanisms for the formation of the seasonal variability of the JES circulation and 

thermohaline structure including the Polar Front meandering and eddies, the TC 

bifurcation and its effect on the formation of mesoscale eddies in the Ulleung Basin and 

the Yamato Basin, the Liman Current (LC) and its penetration into the southwestern 

waters along the Korean coast. The control run, forced by the climatological monthly 
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wind stress and heat flux, is designed to best simulate reality against which each 
experiment is compared. In the experiments, various external and internal fectors are 
modified and the resulting circulation patterns and magnitudes compared to the control 
run results. Specifically, we estimate the contribution, in terms of volume transport and 
circulation patterns, of non-linear advection, wind forcing and lateral boundary transport 
to the ocean features identified in the control results. From this, we can estimate the 
relative importance of these factors to the seasonal variability of the JES circulation and 
thermohaline structure. 
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Fig. 1. The Japan East Sea and surrotmding areas 
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II. JES CURRENT SYSTEMS 


The JES is characterized by a very deep, cold and nearly homogeneous water, the 
Japan Sea Proper Water which occupies the lower part of the basin (Moriyasu, 1972). 
The current coming from the East China Sea enters the JES basin through the 
Tsushima/Korea Straits in the south. The current coming from the Okhotsk Sea enters 
the JES basin through the Tatar Strait in the north; and the current flows out through the 
Tsugaru and Soya Straits in the eastern part of JES (Fig. 2). 

North of 35°N, the Tsushima Warm Current (TC) separates into two branches into 
a western and an eastern channel and flows through the western channel, called the East 
Korea Warm Current (EKWC), and closely follows the Korean coast until it separates 
near into two branches. The eastern branch follows the Polar Front to the western 
coast of Sapporo Island, and the western branch moves northward and forms a cyclonic 
eddy at the Eastern Korean Bay (EKB). It flows through the eastern channel which 
closely follows the Japanese Coast, called the Nearshore Branch (NB) by Yoon (1982a, 
b), or the Japan Coast Current (JCC) by Chu et al. (1999), and it is weaker than through 
the western channel. The strength of the Tsushima Current at both channels reduces with 
depth. 

The cold water enters the JES from the Okhotsk Sea through the Tatar Strait and 
forms the Liman Current (LC) or the North Korean Cold Current (NKCC) carrying 
relatively fresh water along the Russian and North Korean coasts (Seung and Yoon, 
1995). The NKCC meets the EKWC at about 37°N with some seasonal meridional 

migration. After separation from the coast, the NKCC and the EKWC converge to a 
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strong front that runs in a west-east direction across the basin. The NKCC makes a 
cyclonic recirculation gyre in the north but most of the EKWC flows out through the 
outlets (Uda, 1934). The formation of NKCC and separation of EKWC are due to local 
forcing by wind and buoyancy flux (Seung, 1992). Large meanders develop along the 
front and are associated with warm and cold eddies. 

Seung (1995) identified major features of the volume transport from earlier 
numerical modeling results. The transport pattern is largely determined by the upper 
layer circulation and characterized by a large-scale cyclonic recirculation gyre, in which 
the EKWC and the JCC take part, as the inflow-outflow system, and the NKCC in the 
North. At a few hundred kilometers off the separation area, the EKWC makes an 
anticyclonic gyre. The gyre becomes stronger as the EKWC develops. On the other 
hand, the northern cyclonic gyre is very deep and is most significantly in the winter 
strengthened by the wind and buoyancy flux. The gyre, or the southward coastal current 
related to it, is deep enough to intrude southward beneath' the EKWC most of the time. 
Seung also confirmed the summertime presence of counter-current beneath the JCC. 
North of the Polar Front there exists a cyclonic gyre in the Japan Basin (JB) usually 
called the JB gyre. 
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Fig. 2. Surface currents in the JES. EKWC: East Korea Warm Current, PF: Polar Front, 
SWC: Soya Warm Current, TC: Tsushima Current, TWC; Tsugaru Warm Current. From 

Tomczak and Godfrey (1994). 
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in. SEASONAL VARIATION OF THE ATMOSPHERIC SURFACE 

FORCING 


A. GENERAL DESCRIPTION 

The Asian monsoon strongly affects the thermal structure of the JES. During the 
winter monsoon season, a very cold northwest wind blows over the JES (Fig. 3a) as a 
result of the Siberian High Pressure System. By late April, numerous frontally-generated 
events occur making late April and May highly variable in terms of wind speeds and 
number of clouds. During this period storms originating in Mongolia may cause strong, 
warm westerlies (Fig. 3b). By late May and early Jime, the summer surface atmospheric 
low pressure system begins to form over Asia. Initially this low pressure system is 
centered north of the Yellow Sea (YS) producing westerly winds. In late June, this low 
begins to migrate to the west setting up the southwest monsoon that dominates the 
summer months. The winds remain variable through June until Manchurian low pressure 
system strengthens. Despite the very active weather systems, the mean surface wind 
speed over the JES in summer (Fig. 3c) is between 3 and 4 m/s, which is weaker than in 
winter (Fig. 3a). By July, however, high pressure (the Bonin High) to the south and the 
low pressure over Manchuria produce southerly winds carrying warm, moist air over the 
ECS/YS. The simimer monthly mean surface air temperature (SAT) is almost the same 
as the mean sea surface temperature (SST) (Van Loon, 1984). The warm air plus the 
strong downward net radiation stabilizes the upper layer of the water and causes the 
surface mixed layer to shoal. October is the beginning of the transition back to winter 
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conditions. The southerly winds weaken and let the sea surface slope reestablish the 


winter pattern. 

Here, we present a climatological description of the surface net heat and fresh 
water fluxes over the JES. The datasets used were the objectively analyzed fields of 
smface marine climatology and anomalies of fluxes of heat, momentum, and fresh water. 
The fields are derived from individual observations in the Comprehensive Ocean- 
Atmosphere Data Set (COADS) from 1945 to December 1989 and are analyzed on a 1 by 
1 grid (da Silva et al, 1994). The climatological monthly mean values of surface net heat 
and fresh water fluxes are shown in Figs. 4 and 5. 

B. NET SURFACE HEAT FLUX 

Net surface heat flux is computed by 

QNet= Rs-(Rl + Ql + Qs) ( 1 ) 

where Rs is the net downward shortwave radiation, RL the net upward longwave 
radiation, Ql flie sensible heat flux, and Qs the latent heat flux. Positive (negative) values 
of Qwet indicate net heat gain (loss) of the ocean at the surface. The summer field is 
relatively homogeneous (140-160 W m‘^) in the whole JES, whereas a significant 
horizontal gradient increasing from the southeast (Japan coast) to the northwest (east 
Russian coast) exists for the rest of the seasons (Fig. 4). The ocean surface near the 
Tsushima/Korea Strait has the maximum heat loss of 400 W m'^ in the winter (January) 
and the minimum heat gain of 60 W m'^ in the spring (April). This long term net surface 
heat loss will be compensated by the advection of warm waters from the East China Sea. 
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C. SURFACE FRESH WATER FLUX 

The surface fresh water flux is the difference between precipitation rate (P) and 
evaporation rate (E), 

F=P-E (2) 

Positive values of F indicate net water mass gain of the ocean at the surface. The surface 
fresh water flux exhibits a distinct four-season pattern. The winter is featured by fresh 
water gain (2-6 cm/month) in the northern and northeastern JES and fresh water loss (2- 
10 cm/month) in the southern and southwestern JES. A strong horizontal F-gradient is 
found monotonically decreasing from northeast to southwest. The spring (Fig. 5b) and 
summer (Fig. 5c) are both characterized by fresh water gain in the whole JES with 
different horizontal F-gradients: decreasing (increasing) from 4 cm/month (4 cm/month) 
in the northeast JES to 2 cm/month (6 cm/month) in the southwest JES in the spring 
(summer). The autumn (Fig. 5d) is featured by fresh water loss in the whole JES (4-16 
cm/month) with the maximum loss of 16 cm/month near the Tsushima/Korea Strait. 


13 



Latitude (N) 





















































Longitude (E) 


Longitude (E) 


Fig. 5. Precipitation minus evaporation for (a) January, (b) April, (c) July, and (d) 

October, using the COADS ^ta. 

























IV. THE NUMERICAL OCEAN MODEL 


A. MODEL DESCRIPTION 

Coastal oceans and semi-enclosed seas are characterized by extremely high 
temporal and spatial variability that challenges the existing prognostic capabilities of 
numerical simulations. The POM is a three-dimensional, free surface, time dependent, 
primitive equation circulation model on a grid that includes realistic topography 
(Blumberg and Mellor, 1987). Tidal forcing was not included in this application of the 
model, since high frequency variability of the circulation is not evaluated. However, the 
seasonal variations in sea surface height, temperature, salinity, circulation and transport 
are well represented by the model. From a series of numerical experiments, the 
qualitative and quantitative effects of non-linearity, wind forcing and lateral boundary 
transport on the JES are analyzed; yielding considerable insight into the external factors 
affecting the region oceanography. The model results were sampled every thirty days. 

The model contains 94 x 100 x 15 horizontally fixed grid points. The horizontal 
spacing of 10 minutes (approximately 11.54 to 15.18 km in the zonal direction and 18.53 
km in the latitudinal direction) and 15 vertical sigma coordinate levels. 

The model domain is from 35.0° N to 51.5° N, and from 127.0° E to 142.5° E. 
The bottom topography is the smoothed data (Fig. 6) from the Naval Oceanographic 
Office Digital Bathymetry Data Base 5 minute by 5 minute resolution (DBDB5). The 
horizontal diffusivities are modeled using the Smagorinsky (1963) form with the 
coefficient chosen to be 0.2 for this application. The bottom stress Xb is assumed to 
follow a quadratic law 
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U = yOoColVilVi (3) 

where po (= 1025 kg/m^) is the characteristic density of the seawater, V* is the 
horizontal component of the bottom velocity, and Cd is the drag coefficient which is 
specified as 0.0025 (Blumberg and Mellor, 1987) in our model. 

B. ATMOSPHERIC FORCING 


The atmospheric forcing for the JES application of the POM includes mechanical 


and thermohaline forcing. The wind forcing is depicted by 


p oKm 


{du dv^ 
dz ' dz j 


= 0 = (r ojf, T 


( 4 ) 


where {u, v) and (xox, xoy) are the two components of the water velocity and wind stress 
vectors, respectively. The wind stress at each time step is interpolated from monthly 
mean climate wind stress from Comprehensive Ocean-Atmosphere Data Set (COADS), 
1945-1989, with the resolution of 1® x 1°. We interpolated the COADS wind stress to the 
model grid with the resolution of 10’. 

Surface thermal forcing is depicted by 


+ a iC {^OBS — ^ ) (5) 

-aiFS + a^CfSoBs -S) (6) 

where 6bBs and Sobs are the observed potential temperature and salinity, and Cp is the 
specific heat. The relaxation coefficient C is the reciprocal of the restoring time period - 
for a xmit volume of water. The parameters (ai,a 2 ) are (0,l)-type switches: ai = 1, 02 = 
0, would specify only flux forcing is applied; ai = 0, a 2 = 1, would specify that only 


Kh -= a 

dz 


Qhlet 

pCp 
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restoring type forcing is applied. 

In this study, the surface thermohaline forcing is determined solely by the flux 
forcing, that is ai = 1 and a 2 = 0 in (5X6). The mixing coefficients Km,Kh, and Ks were 
computed using a level two-turbulence closure scheme (Mellor and Yamada, 1982). 

C. LATERAL BOUNDARY FORCING 

Closed lateral boundaries, i.e., the modeled ocean bordered by land, were defined 
using a free slip condition for velocity and a zero gradient condition for temperature and 
salinity. No advective or diffusive heat, salt or velocity fluxes occur through these 
boundaries. 

At open boundaries, the numerical grid ends but the fluid motion is unrestricted. 
Uncertainty at open boundaries makes marginal sea modeling difficult. Three 
approaches, local-type, inverse-type, and nested basin/coastal modeling, are available for 
determining the open boundary condition. Here, we take the local-type approach, i.e., to 
use the radiative botmdary condition with specified volume transport. When the water 
flows into the model domain, temperature and salinity at the open botmdaiy are likewise 
prescribed from the climatological data (Levitus, 1982). When water flows out of the 
domain, the radiation condition was applied, 

+ = 0 (7) 

where the subscript n is the direction normal to the boundary. 

Warm water enters the JES through the Tsushima Strait with the Tsushima 
Current from the East China Sea, and cold water enters the JES through the Tatar Strait 
with the Liman Current from the Sea of Okhotsk (Tomczak and Godfrey, 1994). The 
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water exits the JES through the Tsugaru and Soya straits. We found different estimations 
for the flows through three of the JES straits. It was not an easy job to choose one among 
them. There is not much information about the transport through the Tatar Strait, which 
was taken as 0.2 Sv in this study. We use the monthly mean volume transport, reported 
by Yi (1996), through the Tsushima Strait with the annual average of 1.3 Sv, a maximum 
of 2.2 Sv in October, and a minimum of 0.3 Sv in February. The total inflow transport 
through the Tatar and Tsushima straits should be the same as the total outflow transport 
throu^ the Tsugaru and Soya straits. We assume that 75% of the total inflow transport 
should flow out of the JES through the Tsugaru Strait, and 25% throu^ the Soya Strait. 
T his ratio is adopted from the maximum volume transport through the Tsugaru Strait 
estimated by Toba et al. (1982), and through the Soya Strait estimated by Preller and 
Hogan (1998). The monthly volume transports at open boundaries are listed in Table 1. 
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Month 

Jan. 




May 

Jun. 

Tatar Strait(inflow) 

0.2 

Soya Strait (outflow) 

-0.32 

-0.12 

-0.20 

-0.15 

-0.20 

-0.35 

Tsugaru Strait(outflow) 

-0.68 

-0.38 

-0.60 

-0.45 

-0.60 

-1.05 

Tsushima Strait(inflow) 

0.7 

0.3 

0.6 



1.2 


Month 

Jul. 

Aug. 

Sep. 

Oct. 

Nov. 

Dec. 

Tatar Strait(inflow) 

0.2 

Soya Strait (outflow) 

-0.55 

-0.57 




-0.4 

Tsugaru Strait(outflow) 

-1.65 

-1.73 



mi 

-1.2 

Tsushima Strait(inflow) 

2.0 

2.1 

2.1 

2.2 


1.4 


Table 1. The monthly variation of volume transport (1 Sv = 10^ mVs). Here positive 
(negative) values mean inflow (outflow). 


D. INITIAL CONDITIONS AND INITIALIZATION 

The model was integrated with all three components of velocity (u, v, w) initially 
set to zero, and with temperature (Fig. 7) and salinity (Fig. 8) specified by interpolating 
climatological data (Levitus, 1982) to each model grid point. There is no evident Polar 

Front identified from this climatological dataset. The model year consists of 360 days 
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(30 days per month). In order to first capture the winter monsoon, the model was started 
from day 330 (30 November); day 361 corresponds to 1 December. The control run was 
integrated from day 330 to day 1440. The sensitivity runs were integrated from day 
720’s control run results for two years. We use the fourth year’s results to compare the 
JES circulation between the sensitivity and control runs. 

E. MODE SPLITTING 

For computational efficiency, the mode splitting technique (Blumberg and Mellor, 
1987) is applied. The barotropic mode has the time step of 24 sec, based on the Courant- 
Friederichs-Levy (1928) computational stability (CFL) condition and the external wave 
speed. The baroclinic mode has a time step of 720 sec, based on the CFL condition and 
the internal wave speed. 

F. EXPERIMENT DESIGN 

Our approach was to carry out four numerical experiments: one control and three 
sensitivity runs. All runs were completed for the same two-year period encompassing 
both summer and winter monsoons, and, except as specified below, utilized the same 
initial conditions (on day 331). 

Run 1 is the control run. The three sensitivity runs are Run 2, linear dynamics; 
Run 3, no wind; and Run 4, no lateral transport at the open boundaries. The difference 
between the control and the sensitivity runs at each grid point and time should then 
isolate the nonlinear dynamics, lateral boundary effect, and wind forcing effect, allowing 
independent examination of each. 

The theimohaline structure, circulation patterns, and volume transport that 
constitute the JES oceanography will be identified from the control run for subsequent 
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comparative analysis. An important assumption is made that the differences are linear, 
i.e., higher order terms and interactions are negligible and can be ignored. This 
assumption will be examined and qualified. 
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Fig. 6. Smoothed bathymetry 
















































































V. TBDE SIMULATION (CONTROL RUN) 


A. TEMPERATURE 

Fig. 9 shows the sea surface temperature (SST) simulated at day-390 (30 
January), day-480 (30 April), day-570 (30 July), and day-660 (30 October), respectively. 
Comparing to the initial temperature field (Fig. 7), the model successfully simulated the 
formation of the JES Polar Front. Although SST field (Fig. 9) shows an evident seasonal 
variation, the Polar Front exists at all times throughout the year. 

Latitudinal cross-section (135°E) of the temperature field (Figs. 10 and 11) shows 
a strong north-south thermal asymmetry across the Polar Front. As the depth deeper than 
500-m, the temperature is uniformly cold (1-2°C). In the winter (Fig. 10a), the 
thermocline appears in the southern JES, and disappears north of 41°N. The strength 
decreases with latitude from 7"C/100m near the Japan coast to less than Uc/lOOm at 
40°N. From winter to spring (Fig. 10b), the thermocline weakens which however, occurs 
in the whole JES. We still see the northward decreasing in strength. In the summer, a 
strong seasonal thermocline, caused by the strong surface heating, is simulated in the 
shallow depths overlying the permanent thermocline (Fig. 11a). The seasonal 
thermocline sustains from summer to fall (Fig. 11b). We notice that the simulated 
vertical gradient of the seasonal thermocline is too large, which indicates that the vertical 
heat diffusion might not be sufficient in the POM. 
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B. SALINITY 


Three important features are found from the monthly sea surface salinity (SSS) 
field (Fig. 12): (1) salinity is higher in the winter than in the summer; (2) the seasonal 
variability of the surface salinity field is stronger in the south than in the north of the 
Polar Front; and (3) there are two saline activity centers located at the Tatar Strait and the 
Tsushima Strait, respectively. 

Latitudinal cross-sections (135°E) of salinity show a strong north-south haline 
asymmetry across the Polar Front (Figs. 13 and 14): mid-level (200-500 m) salty tongue 
(34.1 psu) appears south of the Polar Front and disappears north of the Polar Front. The 
shallow and strong halocline associated with strong thermocline makes the sub-surface 
water mass north of the Polar Front hydrostatically stable. The mid-level salty tongue 
associated with weak and wide thermocline makes the water mass hydrostatically less 
stable south of the Polar Front. 

C. CIRCULATION 

Fig. 15 shows the sea surface velocity vectors simulated at day-1140 (30 January, 
winter), day-1230 (30 April, spring), day-1320 (30 July, summer), and day-1410 (30 
October, fall), respectively. The simulated surface velocity field (Fig. 15) coincides with 
the earlier description of the JES circulation presented in Chapter 2. The TC separates at 
the Tsushima/Korea Strait into two branches into a western and an eastern channel. Flow 
through the western channel (i.e., EKWC) closely follows the Korean coast until it 
separates near 38°N into two branches. The eastern branch follows the Polar Front to the 
western coast of Sapporo Island, and the western branch, moves northward and forms a 

cyclonic eddy at the Eastern Korean Bay (EKB). Flow through the eastern channel (i.e., 
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JCC) is a little weaker than through the western channel. The LC carries fresh and cold 
water along the Russian coast and becomes the NKCC at the North Korean coast. The 
simulated NKCC meets the EKWC at about 38°. After separation from the coast, the 
NKCC and the EKWC converge to a strong front that runs in the west-east direction 
across the basin. The seasonal variation of the JES circulation is also simulated. 

To investigate the JES circulation as a whole, we integrate the velocity vertically 
from the surface (z = 0) to the bottom (z = -h) 

U=\udz, V==fvdz (8) 

-A -h 


Due to the continuity the volume transport stream function (P) is defined by 


U = - 




( 9 ) 


and satisfies the Poisson equation 


^2 dV dU 

vv =- 

dx dy 


( 10 ) 


For each time instance, we solve the two-dimensional Poisson equation (10) with the 

given boundary conditions, that is, zero transport at the rigid boundaries and known 

transport at the open boundaries (see Table 1). 

Fig. 16 shows the contom plots of the volume transport stream function simulated 

at day-1140 (30 January, winter), day-1230 (30 April, spring), day-1320 (30 July, 

summer), and day-1410 (30 October, fall), respectively. The unit is in Sv (1 Sv = 10^ 
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mVs). The volume transport moves with the smaller values on the left.. The strength of 
the transport is evaluated by the difference between two contours. The northern JES is 
occupied by a cyclonic gyre (called the JB gyre) and the southern JES is characterized as 
a multi-eddy structure. The volume transport streamfimction ('F) has a double-gyre 
structure with negative values (cyclonic) in the northern JES and positive values 
(anticyclonic) in the southern JES. It has an evident seasonal variation with the 
appearance of the cyclonic gyre in the northern JES (i.e., the JB) in the winter and spring 
and disappears in the summer and fall. 

1. Japan Basin (JB) Cyclonic Gyre 

A large-scale cyclonic recirculation gyre over the JB is simulated with a strong 
seasonal variation. This gyre is easily identified by the volume transport streamfimction 
(Fig. 16). The JB gyre is the strongest and recirculates 8 Sv in the winter. It weakens 
and retreats northward in the spring and summer. In the fall, the cyclonic JB gyre 
disappears and weak anticyclonic eddies appears. Our simulation is consistent with the 
earlier study (Seung and Yoon, 1995). 

2. Liman Current (LC) 

The LC is a southwestward current following along the Russian coast. It 
bifurcates into two branches near the east Russian bight (134°E, 42'’N): the western 
branch flows along the Russian-Korean coast and becomes the NKCC. The eastern 
branch flows southeastward, then turns eastward at 41.5°N, and becomes the south frank 
of the Japan Basin (JB) gyre. The LC has a strong seasonal variation with the maximum 
speed in the winter (Fig. 15a) and the minimum speed in the summer (Fig. 15c). 
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Zonal cross-sections of the v-velocity at 46°N for four different seasons indicate 
seasonal and spatial variabilities of the LC. It has a maximum southward component 
(0.21 m/s), occurring near the surface in winter (Fig. 17b) with the width of 100 km and 
the depth of 800 m. The core of the LC is close to the coast and near the surface. In the 
spring, its strength weakens to 0.18 m/s (Fig. 18), but the width and depth keep almost 
unchanged. It further weakens to a minimum of 0.15 m/s in the summer (Fig. 19b) and in 
the fall (Fig. 20b), and shrinks its size into a width of 60 km and depth of400 m. 

3. East Korean Warm Current (EKWC) 

EKWC, western branch of TC north of the Tsushima/Korea Strait, acts as a 
western boimdaiy current and has a strong seasonal variability. At 37®N, the EKWC v- 
velocity component, which is practically the direction of the current on that position, 
varies from 0.42 m/s in summer (Fig. 21a), to 0.30 m/s in winter (Fig. 22a). The width of 
EKWC is around 60 km all year round. However, the depth is around 1,400 m in the 
summer and 800 m in the winter. 

4. EKWC - NKCC Confluence 

The (northward) overshooting of the EKWC near the Korean Bight at 37.5° N is 
stronger in the winter than in the summer (Fig. 15). The overshot EKWC leaves the 
Korean coast and moves northward. It converges with the southward flowing NKCC, at 
39°N, 130°E, and forms a current meandering toward the east along the Polar Front, 
defined as the Polar Front Current (PFC) by Chu et al. (1999). The resultant velocity of 
the PFC is 0.32 m/s and northeastward, in the fall (Figs. 20a and 23a), and 0.21 m/s in the 
winter (Figs. 17a and 24a). 
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5. Japan Coastal Current (JCC) 

JCC, the eastern branch of TC north of the Tsushima/Korea Strait at 132®E, has a 
maYimnm eastward component of 0.24 m/s in the winter (Fig. 25a). A counter current 
beneath it is simulated with a westward component of 0.12 m/s. In the summer, the 
counter current remains almost unchanged, but the JCC falls down to 0.09 m/s (Fig. 26a). 
Seung and Yoon (1995) observed the counter current in the summer. However, we 
simulated the counter current all year round. 

The effect of coastal geometry, such as the Noto Peninsula, on the JCC is also 
simulated. The JCC current turns its direction from the east^^^rd to north-to- 
northeastward along the Noto Peninsula (Fig. 15). The inflow from the Tsushima/Korean 
Strait has a strong seasonal variation with the maximum transport in the fall and the 
miniiTuim transport in the spring (Table 1). The v-component near the tip of the Noto 
Peninsula should have a corresponding variation. Such a topographic steering is clearly 
seen from the zonal cross section of the v-component at 37”N. The northward component 
has a mayimum (mi nimum ) value of 0.21 (0.15) m/s in fall (spring) near 137°E, as shown 
in Fig. 27. Furthermore, the north-to-northeastward JCC near the Noto Peninsula is quite 
narrow and shallow. The width is around 50 km, and the depth is nearly 300 m. Away 
from the narrow JCC, there exists a counter current further off shore and an undercurrent, 
both southward. These currents are weak in the spring (0.06 m/s) and become evident 
(0.09 m/s) in the fall. 

6. Ulleung Basin (UB) Eddy 

The western branch of the warm TC (i.e., EKWC) moves northward along the 

southern part of Korean coast and separates from the coast after approaching the East 
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Korean Bay (EKB). It keeps its northward motion until meeting the NKCC near 40°N, 
meanders southeastward, and forms a warm-core anticyclonic eddy (Fig. 15). The 
simulated anticyclonic UB eddy is strongest in the summer (Fig. 15c). The center is 
located at 38.5°N 130°E. The size of the eddy is around 150 km. The tangential velocity 
is around 0.4 m/s. 
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Fig. 10. Temperature (® C) results from first year on cross-sections at 135° E for (a) 

January and (b) April. 















Fig. 11. Temperature (° C) results from first year on cross-sections at 135® E for (a) July 

and (b) October. 
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Fig. 12. Surface salinity field from the first year results for (a) January, (b) April, (c) July 

and (d) October. 
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Fig. 14. Salinity (psu) results from first year on cross-sections at 135° E for (a) July and 

(b) October. 
















Fig. 15. Surface circulation from control run for (a) January, (b) April, (c) July and (d) 

October. 
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Fig. 16. Voliime transport for (a) January, (b) April, (c) July and (d) October 
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Fig. 17. V component of velocity (m/s) at latitudinal cross sections for January. 
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Fig. 19. V component of velocity (m/s) at latitudinal cross sections for July. 
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Fig. 23. U component of velocity (m/s) at longitudinal cross sections for October 
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Fig. 27. V component of velocity (m/s) at latitudinal cross sections during (a) October and 











VI. DRIVING MECHANISMS 


We analyzed the results of the four experiments to identify the driving 
mechanisms for the JES circulation. 

A. EFFECTS OF NON-LINEARITY (RUN 1 - RUN 2) 

In the first sensitivity study, the nonlinear advection terms were removed from the 
dynamic equations. Otherwise, the same parameters as the control run were used. The 
figures and analysis were performed after computing the difference between control and 
non non-linear terms runs. The differential surface vector velocities (Fig. 28) does not 
show evident general circulation patterns except couple mesoscale eddies, and especially 
the UB eddy, which indicates that the non-linearily does not change the general 
circulation pattern except for the UB anticyclonic eddy. However, our computation 
shows that the non-linearity causes a noticeable change in the volume transport (Fig. 29). 

Fig. 29 shows the plots of the difference in the volume transport streamfunction 
(T) for the control run minus the linear run and represents nonlinear effects. The major 
feature of 'F is the multi-eddy structure with a strong cyclonic eddy with a minimum 
value of -6 Sv near the NKCC and EKWC confluence (131-133 E, 40-42 N). This 
indicates that the non-linear effect is responsible for the eddy generation. Besides, in the 
spring (Fig. 29b) there is a cyclonic gyre with a minimum value of -3 Sv in the JB, 
indicating that the non-linear effect might contribute nearly 60% to the JB cyclonic gyre. 

B. EFFECTS OF WIND (RUN 1 - RUN 3) 

The difference between Run 1 and Run 3 shows the monsoon wind effects. The 

winter wind effect on the surface circulation (Fig. 30a) and the volume transport (Fig. 
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31a) is in contrast to the summer wind forcing (Figs. 30c and 31c). The winter wind 
forcing causes a strong basin-wide JES cyclonic gyre (Fig. 31a) with 8 Sv recirculation 
and the summer wind forcing drives a weak nearly basin-wide JES anticyclonic gyre 
(Fig. 31c). Similar patterns and values are foimd between 'P (Fig. 31) and 'P for the 
control run (Fig. 16) in the northern JES for all the seasons. This indicates that the wind 
forcing is the most important factor (80%) for the generation of the JB cyclonic gyre. 
Different patterns and values are found between 'P (Fig. 31) and 'Pc (control run) (Fig. 
16) in the southern JES in the winter; negative values for T' and positive values for 'Pc. 
This indicates that the wind forcing in the winter weakens the southern JES circulation. 

Compared to the surface velocity vector map for die control run (Fig. 15), the 
surface differential current vectors AV (Fig. 30) clearly show the wind effects on the JES 
circulation: (a) driving the LC in the winter, (b) damping the EKWC especially in the 
winter, (c) generating the UB eddy in all the seasons, and (d) generating eddies along the 
JCC. The wind has almost no effect on the occurrence of the JCC for all seasons. 

The winter monsoon winds blow from the northwest to the southeast over the JES 
surface (Fig. 3) and drives the Ekman flow in the upper ocean to the right of the wind 
direction. Such a surface current moves southward at the western coast of the JES and 
strengthens the LC and weakens the EKWC. The LC strengthens in tbe winter by 30% 
(v-component 0.06 m/s southward) due to wind at 46°N (Fig. 32b). The summer 
monsoon winds blow from the south and southeast to the northwest with much smaller 
speeds. Thus, the summer wind effects on reducing LC and strengthening EKWC are 
quite weak. 
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The winds drive mesoscale eddies especially near UB and along the west coast of 
Honshu. The surface UB eddy, occurring on the differential velocity vector map (Fig. 
30), has the similar swirl velocity (anticyclonic) with same order of magnitudes 
(maximum value of 0.7 m/s) as the control run (Fig. 14). This means that the wind effect 
is in fact a key factor for generating the UB anticyclonic eddy. Furthermore, there is no 
evident current along the west coast of Honshu on the AV map (Fig. 30), indicating that 
wind forcing is not a major factor for maintaining JCC. 

C. EFFECTS OF LATERAL BOUNDARY TRANSPORT (RUN 1 - RUN 4) 

The third sensitivity study used the control run equations and forcing but closed 
all open lateral boundaries, preventing the transport of mass, heat or salinity through the 
Tsushima, Tsugaru, Soya, and Tatar Straits. With inflow or outflow, the JES vertically 
integrated circulation pattern is more pronounced in the summer than in the winter. This 
is because T (control - no transport) pattern (Fig. 33) is quite similar to 'Pc (control run) 
(Fig. 16) in the summer than in the winter. Increased circulation in the summer by the 
lateral transport generally leads to greater horizontal and vertical variabilities of the 
current structure. In the winter, ^ has positive values (2 Sv), and 'Pc has negative values 
in JB. This suggests that the lateral transport reduces the JB cyclonic gyre by 2 Sv in the 
winter. 

Compared to the surface velocity vector map for the control run (Fig. 15), the 
surface differential current vectors AV (Fig. 34) clearly show the lateral transport effects 
on the JES circulation: (a) driving the UB eddy in all the seasons, (b) generating EKWC, 
and (c) generating JCC and eddies along the JCC. 
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An eddy over the Ulleung Basin (38°N 130°E) is well depicted in the summer 
surface AV field (Fig. 34c). The calculated swirl velocity is 0.25 m/s at 39°N 130®E 
(Figs. 35a and 36a); similar magnitude to what was observed in the non-linear effect. 
However, the depth ranges from 300-600 m, vvhich are shallower than the non-linear 
case. 

Along the East Korean Coast in the winter, Av has a value of 0.15 m/s at 37®N 
(Fig. 37a), 50% of the v-component in the control run. Finally, at 40°N 130°E the 
calculated |AV1 is 0.21 m/s to the northeast direction (Figs. 37b and 38a), reaching a 
depth of 1400m. This corresponds to 100% of |V| in the control run (Figs. 22b and 24A). 
This suggests that in the winter the lateral transport contributes 50% (100%) to the 
EKWC in the southern (northern) Korean Coast. In the summer, Av has a velocity of 
0.09 m/s at 37°N (Fig. 39a), or 21% of the v-component in the control run. 

The lateral transport has very weak effects on the occurrence of the LC current for 
all seasons. The LC strength in the winter is 30% due to the lateral transport at 46°N 
(Fig. 40b). 
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Fig. 29. Volume transport from non-linear effects for (a) January, (b) April, (c) July and 

(d) October. 
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Fig. 30. Surface circulation from no wind effects for (a) January, (b) April, (c) July and 

(d) October. 


















































































































































Longitude (E) (at Latitude (N) = 46) 


Fig. 32. V component of velocity (m/s) at latitudinal cross sections from no wind effects 

for January. 
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Fig. 33 Volume transport from no lateral transport effects for (a) January, (b) April, (c) 

July and (d) October. 
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Fig. 34. Surface circulation from no lateral transport effects for (a) January, (b) April, (c) 

July and (d) October. 































































































Longitude (E) (at Latitude (N) = 46) 


Fig. 35. V component of velocity (m/s) at latitudinal cross sections from no lateral 

transport effects for July. 





















Longitude (E) (at Latitude (N) = 40) 



Fig. 37. V component of velocity (m/s) at latitudinal cross sections from no lateral 

transport effects for January. 
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Fig. 38. U component of velocity (m/s) at longitudinal cross sections from no lateral 

transport effects for January. 
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Fig. 39. V component of velocity (m/s) at latitudinal cross sections from no lateral 

transport effects for July. 
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VII. CONCLUSIONS 


(1) The JES circulation and thermohaline structure was simulated in this study by 
the POM model under the climatological forcing. The northern JES is occupied by a 
cyclonic gyre (called the JB gyre) and the southern JES is characterized as a multi-eddy 
structure. The volume transport streamfunction has a double-gyre structure with negative 
values (cyclonic) in the northern JES and positive values (anticyclonic) in the southern 
JES. It has an evident seasonal variation with the appearance of the cyclonic gyre in the 
northern JES (Le., the JB) in the winter and spring and disappears in the summer and fall. 

(2) A large-scale cyclonic recirculation gyre over the JB is simulated with a 
strong seasonal variation. The JB cyclonic gyre is the strongest and recirculates 8 Sv in 
the winter. It weakens and retreats northward in the spring and summer. In the fall, the 
cyclonic JB gyre disappears and weak anticyclonic eddies appear. 

(3) The POM model simulates the basic currents in the JES such as the LC, 
EKWC, NKCC, and JCC. The LC has a maximum southward component (0.21 m/s) and 
occurs near the surface in the winter with a width of 100 km and a depth of 800 m. The 
core of the LC is close to the coast and near the surface. In the spring, its strength 
weakens to 0.18 m/s, but the width and depth remain almost unchanged. It further 
weakens to a minimum of 0.15 m/s in the summer and fall, and shrinks in size to a width 
of 60 km and depth of 400 m. The EKWC, the western branch of the TC north of the 
Tsushima/Korea Strait, acts as a western boundary current varying from 0.42 m/s 
(summer) to 0.30 m/s (winter). The width of the EKWC is around 60 km all year round. 
However, the depth is around 1,400 m in the summer and 8(K) m in the winter. The 
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(northward) overshot EKWC leaves the Korean coast and moves northward. It converges 
with the southward flowing NKCC, at 40°N, 130°E, and forms a current meandering 
toward east along the Polar Front, defined as the Polar Front Current (PEC). The 
resultant velocity of the PEC is 0.32 m/s and northeastward, in the fall, and 0.21 m/s in 
the winter. 

The JCC, the eastern branch of the TC north of the Tsushima/Korea Strait at 
132°E, has a maximum eastward component of 0.24 m/s in the winter. A counter current 
beneath it is simulated with the westward component of 0.12 m/s. In the summer, the 
counter current keeps almost unchanged, but the JCC falls down to 0.09 m/s. This is 
consistent with the observational study reported by Seung and Yoon (1995). The effect 
of coastal geometry, such as the Noto Peninsula on the JCC, is also simulated. 

(4) The UB warm-core anticyclonic eddy is also simulated. The western branch 
of the warm TC (ie., EKWC) moves northward along the southern part of the Korean 
coast and separates from the coast after approaching the East Korean Bay (EKB). It 
keeps its northward motion until meeting the NKCC near 40°N, meanders southeastward, 
and forms a warm-core anticyclonic eddy. The simulated anticyclonic UB eddy is 
strongest in the summer. The center is located at 38.5°N 130°E. The size of the eddy is 
around 150 km. The tangential velocity is around 0.4 m/s. 

(5) The POM simulates the formation of the JES Polar Front and its seasonal 
variation. The simulated thermal field shows a strong north-south thermal asymmetry 
across the Polar Front. At a depth deeper than 500 m, the temperature is uniformly cold 
(1-2°C). In the winter, the thermocline appears in the southern JES, and disappears north 
of 41®N. The strength decreases with latitude from 7°C/100m near the Japan coast to less 


72 




than l°C/100m at 40°N. In the summer, a strong seasonal thermocline, caused by strong 
surface heating, is simulated in the shallow depths overlying the permanent thermocline. 
The seasonal thermocline is sustained from summer to fall. 

(6) The POM simulates the salinity field reasonably well with the surface featured 
by (a) higher salinity in the winter than in the summer; (b) stronger seasonal variability in 
the south than in the north of the Polar Front; and (c) two salinity activity centers located 
at the Tatar Strait and the Tsushima Strait. The model also simulates a strong north-south 
salinity asymmetry across the Polar Front with the salty tongue (34.1 psu) appearance 
only in the south of the Polar Front. The shallow and strong halocline associated with 
strong thermocline makes the sub-surface water mass north of the Polar Front 
hydrostatically stable. The mid-level salty tongue associated with weak and wide 
thermocline makes the water mass hydrostatically less stable south of the Polar Front. 

(7) The nonlinear advection does not affect the general circulation pattern 
evidently, but does affect the formation of the mesoscale eddies, and especially the UB 
eddy (aU seasons) and the JB gyre (spring). 

(8) The model wind effects on the JES circulation are more pervasive than those 
of non-linear dynamic effects. The winter winds cause a strong basin-wide JES cyclonic 
gyre with 8 Sv recirculation in the northern JES and the summer winds drive a weak 
nearly basin-wide JES anticyclonic gyre. Thus, the wind forcing is the most important 
fact (80%) for the generation of the JB cyclonic gyre. The winds also influence the 
surface circulation such as driving the LC (winter), damping the EKWC (winter), 

generating the UB eddy (all seasons), and generating eddies along the JCC (aU seasons). 
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The wind has almost no effect on the occurrence of the JCC for all seasons. 

(9) The model boundary-forcing enhances (weakens) the JES volume transport in 
the summer (winter). It has very weak effects on the occurrence of the LC except in the 
winter, when the boundary-forcing accounts for 30% of the LC at 46°N. It weakens the 
JB cyclonic gyre by 2 Sv (25%) in the winter. Besides, the boundary forcing also 
influences the surface circulation such as driving the UB eddy in all the seasons (50% in 
the winter), generating the EKWC (50-100% in the winter and 21% in the summer), and 
generating the JGC and eddies along the JCC. 

(10) Future studies should concentrate on less simplistic scenarios. Realistic 
lateral transport should be included and the use of extrapolated climatological winds 
needs to be upgraded to incorporate synoptic winds to improve realism. Finally, the 
assumption of quasi-linearity that allowed us to use simple differences to quantify the 
effect of external forcing needs to be rigorously tested. It is important to develop a 
thorough methodology to perform sensitivity studies under the highly non-linear 
conditions that may exist in the httoral environment. 
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